Dynamically stable multiply quantized vortices in dilute Bose-Einstein condensates 
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Multiquantum vortices in dilute atomic Bose-Einstein condensates confined in long cigar-shaped 
traps are known to be both energetically and dynamically unstable. They tend to split into single- 
quantum vortices even in the ultralow temperature limit with vanishingly weak dissipation, which 
has also been confirmed in the recent experiments [Y. Shin et al, Phys. Rev. Lett. 93, 160406 
(2004)] utilizing the so-called topological phase engineering method to create multiquantum vortices. 
We study the stability properties of multiquantum vortices in different trap geometries by solving the 
Bogoliubov excitation spectra for such states. We find that there are regions in the trap asymmetry 
and condensate interaction strength plane in which the splitting instability of multiquantum vortices 
is suppressed, and hence they are dynamically stable. For example, the doubly quantized vortex 
can be made dynamically stable even in spherical traps within a wide range of interaction strength 
values. We expect that this suppression of vortex-splitting instability can be experimentally verified. 



I. INTRODUCTION 

Quantized vortices are fundamental topological exci- 
tations in systems having long range quantum phase co- 
herence, e.g., in superfluids and superconductors. Cre- 
ation and observation of singly quantized vortices in di- 
lute atomic Bose-Einstein condensates (BECs) in 1999 [J 
was an important demonstration of the superfluid prop- 
erties of these gaseous systems. Since then, the proper- 
ties of singly quantized vortices and vortex lattices have 
been studied extensively, both experimentally and theo- 
retically 0. 

However, there has been considerably fewer investiga- 
tions of multiply quantized vortices, for which the phase 
of the order parameter winds a multiple k of 27r. Typ- 
ically, quantized vortices are generated by rotating the 
condensate with a time-dependent external potential. 
Since the energy of multiquantum vortices in harmoni- 
cally trapped condensates increases considerably faster 
than linearly with respect to vorticity quantum number 
K, single-quantum vortices are favourable when the sys- 
tem is rotated and thus multiquantum vortices do not 
appear. 

Although difficult to create with the usual rotational 
techniques, Leanhardt et al. were able to create the first 
two and four times quantized vortices in dilute BECs 
using the so-called topological phase engineering method 
proposed by Nakahara et al. 0,IEIEI3|- In this technique, 
the bias field of the loffe-Pritchard trap is reversed almost 
adiabatically, which transfers the atoms coherently from 
the original irrotational spin state to other spin states 
that contain multiquantum vortices. In the vicinity of 
the middle point of the field reversal period the bias field 
is weak, and hence spin flips may occur close to the center 
of the trap. Due to these Majorana flips approximately 
half of the atoms are transferred into nonconflned states 
and escape from the trap. If the total spin quantum 
number of the condensate atoms is = 1, the bias field 
reversal transfers atoms into the single weak field seeking 
state (WFSS) which contains a doubly quantized vortex. 



On the contrary, if the condensate consists of atoms with 
F = 2, there are two WFSSs. In the theoretical sim- 
ulations 0, both of these WFSSs are macroscopically 
populated due to Majorana fiips after the field reversal. 
One of the WFSSs contains a three times quantized vor- 
tex and the other a four times quantized one. In the 
experiments however, the measured spatial angular mo- 
mentum per particle was (4.4 ± 0.4)7i, seeming to match 
with the angular momentum of a pure four times quan- 
tized vortex Ij] . Since no tomographic images of different 
spin components for F = 2 condensates after field rever- 
sal has been reported as for the F = 1 condensates Q, 
it still remains somewhat open whether the created state 
was actually a mixture of three- and four-quantum vor- 
tices. 

Theoretical investigations 0, 1^ have shown that mul- 
tiply quantized vortices are dynamically unstable in the 
experimentally realized cigar-shaped geometries 0, 0, 
[ij. This implies that they do not only decay at fi- 
nite temperatures in the presence of dissipation, but also 
small perturbations to a stationary multiquantum vortex 
in a pure condensate result in its splitting into k singly 
quantized vortices. The splitting instability of doubly 
quantized vortices has also been observed in the exper- 
iments pTl |. in which their lifetimes were measured as 
functions of the particle density. It has also been the- 
oretically confirmed that the gravitational sag and the 
time-dependent magnetic potential used in the experi- 
ments alone suffice to act as the perturbation required 
to initiate the splitting process and lead to agreement 
with the measured lifetimes even when dissipation is ne- 
glected JXy. Furthermore, compensation of the gravita- 
tional sag with an optical potential has resulted in longer 
lifetimes of condensates with four-quantum vortices |l2l| . 

To increase the lifetime of multiquantum vortices, one 
could try to minimize perturbations breaking the rota- 
tional symmetry of the vortex state, for example by ar- 
ranging the gravitation to be parallel to the symmetry 
axis of the vortex, or try to choose the system param- 
eters such that the state possibly becomes dynamically 
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stable. Suppression of the rotational perturbations is a 
difScult way to increase the lifetime significantly, since 
the amplitudes of the splitting instability modes begin 
to grow exponentially in time. In the latter possibility, 
the vortex state would not decay in the temporal evolu- 
tion under small perturbations in time scales for which 
dissipative effects can be neglected. 

In this paper, we show that two-, three-, and four- 
quantum vortices in harmonically trapped condensates 
can be made fully dynamically stable by tuning the in- 
teraction strength of the condensate and the aspect ra- 
tio of the atom trap suitably. It is known that in the 
limit of effectively two-dimensional pancake-shaped con- 
densates there exist dynamical stabilit y reg ions with re- 
spect to the interaction strength Il4| . Our results 
show that these stability regions extend nontrivially to 
genuine three dimensional geometries, and for example 
the doubly quantized vortex has wide stability regions as 
a function of the interaction strength even for spherical 
and slightly prolate condensates. The stability regions 
of three- and four-quantum vortices are smaller, but we 
show that also their splitting instabilities can be sup- 
pressed for oblate and spherical condensates. Further- 
more, since the interaction strength parameter can be 
easily tuned by changing the condensate particle number 
or by utilizing atomic Feshbach resonances, we expect 
that this stabilization of multiquantum vortices can be 
verified experimentally. 



II. METHODS 

At ultralow temperatures, the complex valued order 
parameter field ip{r) of a weakly interacting BEG is de- 
scribed by the Gross-Pitaevskii (GP) equation 



-5|V(r,t)|21^(r,t), 



(1) 



where Ti — — + Vi,xt(r), m is the mass of the con- 

fined atoms and T4xt(r) the external potential. The in- 
teraction strength parameter g is related to the vacuum 
s-wave scattering length hy g = 47r ?i^a/m. Throughout 
this paper we consider condensates confined in rotation- 
ally symmetric harmonic traps of the form 



^ 2,' 2 I \2 2\ 

-■muj^(r + A z ), 



(2) 



where r'^ — x'^ + and A = ujzl^r is the ratio between 
the trapping frequencies in the axial and radial direc- 
tions. The order parameter is normalized according to 
/ |'0(r)|^dr = iV, where N is the total particle number. 

Stationary states of the system with an eigenvalue /i 
are of the form ■0(r)e~*''*/''. Small-amplitude oscillations 
about a given stationary state ip{r)e~^f^^^^ can be studied 
by writing the order parameter as 



By substituting this trial function into the GP equation 
with the decomposition 



i?(r,<)-^[u,(r)e— '*+<(r)e'' 
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and linearizing with respect to '!9(r, i), we obtain the Bo- 
goliubov equations 



C{r) M{r) 
-M*{r) -£{r) 



huj„ 



Ug{r) 
Vqir) 



(5) 



where £(r) = H - fi + 2g\ip{r)\'^, M{r) = gtp'^{r), Ug{r) 
and Vq{r) are the quasiparticle amplitudes and ujq is the 
eigenfrequency of the mode q. From Eqs. Q and one 
notes that the existence of non-real eigenfrequencies Uq 
results in exponential initial growth of the population of 
the corresponding mode, i.e., the corresponding station- 
ary state is dynamically unstable. Naturally, in the pres- 
ence of such dynamically destabilizing modes, Eqs. Q 
and Q hold only for the initial dynamics of the conden- 
sate for which the second and higher powers of i9(r, t) are 
negligible. 

Since we consider rotationally symmetric external po- 
tentials given by Eq. the minimal-energy stationary 
state containing a K-quantum vortex can be written as 



(6) 



where (r, 0, z) are the cylindrical coordinates. For this 
kind of stationary state the azimuthal angle dependence 
can be separated from the solutions of the Bogoliubov 
equations lO, and hence the quasiparticle amplitudes can 
be expressed as 



Uq{r) = ^/,(r,z)e^('''+'')^ 



(7) 
(8) 



(3) 



where Kq denotes the angular momentum quantum num- 
ber of the mode. The dynamical stability of a given sta- 
tionary K-quantum vortex state can thus be determined 
by solving effectively two-dimensional Bogoliubov equa- 
tions for all quantum numbers Kq. In practice, however, 
there exist complex eigenfrequencies only for values of Kq 
in the vicinity of k iljj . 

By scaling time with l/tUr, length with a,- = y^h/mujr 
and energy with hujr, one obtains a dimensionless GP 
equation. From this equation, it is observed that the 
stationary K-quantum vortex states are completely de- 
termined by the trap asymmetry parameter A and the 
effective interaction strength g = AnaN/ar. Thus we will 
study the dynamical stability of vortex states as functions 
of these parameters. 



III. RESULTS 

We have solved the stationary states and the corre- 
sponding Bogoliubov excitation spectra for two-, three-, 
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and four-quantum vortices in the parameter ranges < 
1/A < 2 and < ^ < 1000. For typical parameter val- 
ues used in the experiments 0, the value g = 1000 of 
the interaction strength corresponds to the particle num- 
ber iV « 4 • 10"*. To investigate condensates with higher 
particle numbers, we have also calculated the excitation 
spectra for A = 1 and < .g < 5000. 



A. Two-quantum vortex 

Figure ^ shows the maximum imaginary parts of the 
eigenfrequencies uoq for two-quantum vortex states. In 
this case, there is only one angular momentum quan- 
tum number, namely Kq — 2, for which complex frequen- 
cies exist. In the limit 1/A corresponding to thin 
pancake shaped condensates, the imaginary part Im(a;q) 
attains a quasi-periodic form as a function of the inter- 
action strength g. In this limit, the results approach 
those obtained from the effectively one-dimensional cal- 
culations reported in Ref . . As 1/ A increases towards 
unity corresponding to the spherical trap shape, the in- 
stability regions widen with respect to the interaction 
parameter. However, we note that for the two-quantum 
vortex state there are fairly large regions in the param- 
eter space where the vortex state is dynamically stable. 
Only for slightly prolate trap geometries, 1/A « 2, the 
instability regions start to overlap, and the vortex state 
is dynamically unstable for most interaction strengths g 
under consideration. In Fig. 12 the maximum imaginary 
part of the eigenmode frequencies is shown for interac- 
tion strength values in the range < g < 5000, in the 
special case of the spherically symmetric trap with A = 1. 
The largest dynamical stability regions are seen to exist 
at fairly weak interaction strengths. 



B. Three-quantum vortex 

Similarly to the two-quantum vortex case, the maxi- 
mum imaginary parts of the Bogoliubov eigenfrequencies 
are shown in Fig. |31 for stationary three-quantum vortex 
states. In this case, there are three angular momentum 
quantum numbers, Kq = 2,3, and 4, which yield complex 
eigenfrequencies. There exist stable regions also for the 
three-quantum vortex state, but they are considerably 
smaller compared with the two-quantum vortex. The 
large regions of instability are due to modes with quan- 
tum number Kg = 2, while the thinner line-like regions 
originate from modes with Kq = 3. There exist complex 
eigenfrequencies also for Kg = 4, but their absolute value 
is significantly smaller compared with the ones for Kq = 2 
and 3. The similar shape of the instability regions corre- 
sponding to the quantum number Kq — 2 for the two- and 
three-quantum vortex states is also noticeable. Figure 0] 
shows the maximum imaginary parts of the eigenfrequen- 
cies for higher interaction strengths up to g = 5000 for 
three-quantum vortex states in the spherically symmetric 
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FIG. 1: Gray scale plot of the maximum imaginary part of 
the Bogohubov eigenvalues as a function of the interaction 
strength g = AnaN/ar and inverse trap anisotropy 1/A = 
Ur/uJz for the two-quantum vortex state. White denotes the 
dynamically stable regions with vanishing imaginary parts, 
and black corresponds to maximum imaginary part O.lSbur, 
the scale in between being linear. The strip-like instability 
regions for small 1/A are in fact continuous, although they 
seem fragmentary due to discretization of the parameter space 
in the computations. 



0.15 




FIG. 2: Maximum imaginary part of the Bogoliubov eigen- 
values for the two-quantum vortex state as a function of the 
interaction strength g — A-KaN/a,- for the spherically symmet- 
ric trap geometry with A = 1. 



trap geometry. 



C. Four-quantum vortex 

Figure [3 shows the maximum imaginary parts of the 
Bogoliubov eigenfrequencies for the four-quantum vortex 
state. In this case, complex eigenfrequencies exist for ex- 
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FIG. 3: Maximum imaginary part of the Bogoliubov eigen- 
values as a function of the interaction strength g = iiraN/ar 
and inverse trap anisotropy 1/A for the three-quantum vor- 
tex state. In the grayscale white corresponds to dynamically 
stable state and black the imaginary part value 0.208LJr of 
the frequency. The large instability regions are due to modes 
with Kq = 2, and the thin lines in between them are caused 
by modes with Kg — 3. 
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FIG. 4: Maximum imaginary part of the Bogoliubov eigen- 
values for the three-quantum vortex state as a function of the 
interaction strength g = AiraN/ar in the spherically symmet- 
ric trap geometry. 



citation quantum numbers Kq = 2,3,4, and 5. As for 
the two- and three-quantum vortex states, modes corre- 
sponding to the excitation quantum number Kjq — 2 are 
again dominating, i.e., have the largest imaginary parts, 
for most of the parameter values. However, in the limit 
of weak interaction strength, g — > 0, a mode with Kg = 3 
dominates the dynamical instability of the vortex state. 
Figure El shows the maximum imaginary parts for A = 1 
and < g < 5000. There are only a few interaction 
strength values for which the four-quantum vortex is dy- 



namically stable in this geometry. 

Comparison of Figs. 01 and |B1 shows that the dynami- 
cal instability is suppressed for both the three and four 
quantum vortex states for g « 3100. Actually, they are 
both dynamically stable for g = 3135. This stability win- 
dow can possibly be utilized to stabilize the coexisting 
three- and four-quantum vortices created in the topolog- 
ical phase imprinting method for F = 2 condensates. 




FIG. 5: Maximum imaginary part of the Bogoliubov eigen- 
values as a function of the interaction strength g = ^naN/ar 
and inverse trap anisotropy 1/A for the four-quantum vortex 
state. Complex frequencies exist in this case for excitation 
quantum numbers Kq = 2, 3, 4, and 5. In the linear grayscale 
black corresponds to imaginary part value 0.231tJr of the fre- 
quency and the white to zero. 
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FIG. 6: Maximum imaginary part of the Bogoliubov eigen- 
values for the four-quantum vortex state as a function of the 
interaction strength g = 'i-naN/ar in the spherically symmet- 
ric trap geometry. 
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IV. CONCLUSIONS 

We have studied the dynamical stabihty of two-, three-, 
and four-quantum vortex states in dihitc BECs confined 
in rotationally symmetric traps. Experimentally, multi- 
quantum vortex states have primarily been realized in 
long cigar-shaped traps, in which they arc known to be 
dynamically unstable and to split rapidly into single- 
quantum vortices. By calculating the Bogoliubov exci- 
tation spectra of vortex states as functions of the trap 
asymmetry parameter and the condensate particle num- 
ber, we have shown that there exists regions in which 
the splitting instabilities are suppressed and the mul- 
tiquantum vortices are dynamically stable. Especially, 
the stability regions of the two-quantum vortex states 
are large, and extend from oblate pancake-shaped con- 
densates beyond the spherical geometry. For increasing 



vorticity quantum number the number of destabilizing 
modes grows and the stability regions shrink consider- 
ably, but even the four-quantum vortex state can be 
made dynamically stable in the spherical geometry for 
certain condensate particle numbers. We assume that the 
predicted suppression of multiquantum vortex splitting 
instabilities for certain values of trap asymmetry param- 
eter and particle number can be experimentally verified. 
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